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1 .  INTRODUCTION 


A  new  technique  using  multiple  flash  x-ray  sources  for  the  determination 
of  cross-sectional  density  profiles  of  rapidly  translating  or  dynamically 
deforming  objects  Is  presented.  The  method,  an  extension  of  the  concepts  of 
medical  tomography,  allows  stop-action  survey  of  events  with  time  constants  at 
the  low  end  of  the  microsecond  range,  making  It  suitable  for  problems  of 
ballistic  research  and  materials  proof  testing. 

The  rapidity  of  ballistic  events  places  severe  requirements  on  the  design 
of  a  computed  tomography  facility.  In  addition  to  the  short  time  available 
for  data  registration,  spatial  resolution  In  the  millimetre  range  and  a 
density  resolution  better  than  one  percent  Is  required.  Also,  the  x-ray  beam 
has  to  be  energetic  enough  to  traverse  at  least  two  centimetres  of  a  high 
atomic  number  material  such  as  steel  with  the  emerging  beam  having  a  signal  to 
noise  ratio  favorable  for  good  Image  quality.  Such  energy  and  dosage  levels 
are  obtainable  from  commercial  flash  x-ray  units  rated  In  the  megavolt  range. 

In  fourth-generation  conventional  CT  systems,  the  source  Is  rotated 
around  a  stationary  object  with  the  transmitted  radiation  recorded  by  a 
stationary  ring  of  detectors.  Data  acquisition  times  of  a  few  seconds  are 
common.  With  time  constants  of  the  observed  scene  of  the  order  of 
microseconds,  mechanical  motion  of  the  source  Is  no  longer  feasible  and 
another  means  of  obtaining  the  transmission  profiles  must  be  considered. 

The  first  practical  approach  to  CT  of  moving  objects  resulted  In  the 
building  of  the  Dynamic  Spatial  Reconstructor  of  the  Hayo  Clinic.  Ry  means  of 
gating,  and  using  28  separate  thermionic  x-ray  sources,  pictures  of  acceptable 
quality  of  the  beating  heart  have  been  obtained.^  In  view  of  the  long  pulse 
width  (around  300  microseconds)  of  the  thermionic  x-ray  sources,  this  method 
Is  not  suitable  for  ballistic  applications. 

Scanning  electron  beam  systems,  such  as  the  Imatron  0-100  clne-CT 
scanner,^  where  an  electron  beam  Is  magnetically  deflected  along  a  circular 
anode  ring  generating  an  electronically  movable  x-ray  source,  show  promise  of 
higher  temporal  resolution  but  as  yet,  due  to  problems  with  heating  of  the 
anode,  are  limited  to  lower  energies  and  longer  time  scales  than  those 
required  for  ballistic  applications. 

A  solution  to  the  problem  of  generating  all  the  transmission  profiles 
needed  In  a  few  microseconds  Is  possible  by  the  use  of  simultaneously- 
triggered  multi-source  arrangement.  Ideas  along  these  lines  were  advanced  by 


^J.H.  Kinsey,  R.A.  Robb,  E.L.  Rltman,  E.H.  Wood,  "13)0  DSR-  a  High  Temporal 
Resolution  Volumetric  Roentgenographlc  CT  Scanner,"  Hers.  Vol.5,  #3,  pp.  177- 
188,  1980. 

^D.W.  Farmer,  M.J.  Llpton,  C.B.  Higgins,  "Clne-CT  Captures  the  Beating  Heart," 
Diagnostic  Imaging.  Vol.6,  pp.  54-58,  1984. 


Triable  and  Aseltlne^  who  developed  a  six  source  sequentially^fired 
radiographic  facility. 

A  recently  completed  feasibility  study  by  two  of  the  present  authors^  has 
shown  that  industrial  quality  tomographic  data  of  events  having  time  constants 
in  the  microsecond  range  is  indeed  possible.  There,  it  was  demonstrated  that 
on  mock-ups  of  ballistic  interest,  typically  objects  up  to  SO  cm  in  diameter 
and  having  constituents  consisting  of  Si,  C,  N,  0,  a  1  HT  source  would  provide 
transmitted  photons  adequate  in  number  and  energy  to  assure  a  favorable 
signal-  to-noise  ratio  at  the  detector  and  thus  an  acceptable  tomographic 
reconstruction.  Here,  we  discuss  and  expand  on  these  findings  and  detail  the 
design,  the  components,  the  capabilities,  as  well  as  some  of  the  system 
studies  which  led  up  to  the  configuration  which  is  under  construction  and 
scheduled  for  completion  in  the  near  future. 

2.  LAYOUT  OF  THE  FACILITY 

A  schematic  diagram  of  the  tomographic  layout  is  given  in  Figure  1.  The 
twenty-one  sources  and  detector  arrays  (lAich  can  be  either  medical  film- 
screen  cassettes  or  solid  state  detectors)  are  arranged  in  a  360"  arc  about 
the  object.  This  configuration  has  the  advantage  over  the  previous  design^  of 
being  more  compact,  allowing  for  both  more  sources  and  a  smaller  source-to- 
object  distance.  As  a  consequence,  a  larger  dose  will  be  received  by  the 
object  and  the  detector,  increasing  the  sensitivity.  There  are,  however, 
disadvantages.  Shielding  of  the  detectors  from  unwanted  scattered  and  punch- 
through  radiation  will  be  slightly  more  difficult.  Details  of  the  layout  are 
given  in  Table  1. 

Table  1.  System  Parameters 


Scan  diameter  25  cm 

Beam  thickness  2-10  mm 

Number  of  sources  21 

Source  to  object  distance  2.10  m 

Focal  spot  size  1-3  mm 

Pulse  width  25  ns 

Source  voltage  1  W 

Dose  65  mR  0  1  m 

Source  to  detector  distance  3.78  m 

Magnification  variable 


The  system  can  also  be  used  to  obtain  conventional  radiographs.  The 
collimator,  which  determines  the  slice  thickness  (2  to  10  mm),  can  be  removed 
or  widened  to  obtain  21  radiographs.  The  sources  could  be  fired  sequentially 


^J.  Trimble,  C.L.  Aseltine,  "Flash  X-Ray  and  Cineradiography  at  100  000  Fps," 
May  1983,  ARBRL-TR-0249 1 ,  Aberdeen  Proving  Ground,  M),  ADA  129090. 

^C.K.  Zoltanl,  K.J.  White,  R.P.  Kruger,  "Result  of  Feasibility  Study  on 
Computer  Assisted  Tomography  for  Ballistic  Applications,"  August  1983,  ARBRL- 
TR*02513,  Aberdeen  Proving  Ground,  M),  ADA  133214. 


In  tine  for  clne-radlographlc  data.  Additionally,  they  could  be  fired  in 
pairs  giving  10  sets  of  nearly  orthogonal  (85.7")  radiographs.  Sources  could 
also  be  fired  in  adjacent  pairs  to  yield  stereographic  images.  A  potential 
problem  exists  with  this  setup.  The  focal  spot  size  of  conventional  high- 
voltage  flash  tubes  is  5mm.  The  arrangement  shown  in  Figure  1  has  a 
relatively  large  object-to-screen/source-to-object  ratio.  This  will  produce 
radiographs  with  large  penumbras  with  potentially  blurred  Images.  This  can  be 
Improved  by  using  a  demountable  anode  tube  where  smaller  focal  spots  can  be 
realized.  An  alternative  solution  is  to  decrease  the  object-to  screen 
distance.  The  detector  screens  are  moved  onto  a  smaller  circle,  in  Figure 
1.  In  this  position  only  five  radiographs  will  be  produced  by  firing  x-ray 
units  marked  No.  1  through  No.  5.  Alternative  configurations  are  also 
possible  that  will  give  more  than  5  radiographs  but  with  some  sacrifice  of 
resolution. 

The  layout  discussed  up  to  now  allows  the  smallest  object  to  source  and 
object  to  detector  distances.  We  now  discuss  the  advantages  of  this  layout, 
hereafter  referred  to  as  System  III,  see  Figure  1  in  comparison  to  a  previous 
and  also  viable  approach,  described  in  Reference  4,  called  System  II  and  which 
is  shown  in  Figure  2. 

3.  SYSTEMS  ANALYSIS 

3.1  Background 

A  complete  derivation  of  the  equations  which  describe  the  most  compact 
geometries  for  Systems  II  and  III  is  given  in  the  Appendix  A.  These 
geometrical  equations  are  plotted  in  Figure  3  for  a  25  cm  scan  circle 
diameter,  dg  is  the  source-to-object  center  distance  and  d^  is  the  detector- 
to-object  center  distance. 

3 .2  Comparison  of  Systems  II  and  Ill 

It  can  be  seen  directly  that  System  III  is  always  more  compact  than 
System  II.  For  example,  with  a  magnification  of  1.5  and  21  beams,  the  source- 
to-detector  distance  for  System  II  is  7.5  m  whereas  for  System  III  it  is  only 
5.0  m.  The  smaller  size  of  System  III  allows  the  substructure  supporting  the 
sources  and  detectors  to  be  smaller  and  less  costly.  It  further  allows  a 
higher  x-ray  Intensity  at  the  detector  giving  lower  image  noise  and  better 
contrast  resolution. 

Because  in  System  III  the  sources  cover  a  full  360  degrees,  a  more 
uniform  sampling  and  a  balanced  slice  thickness  will  be  obtained.  In  fact. 
System  II,  which  is  constrained  to  have  a  source  arc  of  less  than  180  degrees 
cannot  yield  a  mathematically  complete  reconstruction.  A  complete 
reconstruction  for  a  fan  beam  geometry  requires  a  source  arc  of  180  degrees 
plus  the  fan  angle. 


The  final  advantage  of  System  III  relates  to  the  fact  that  it  has  nearly 
opposing  detectors.  By  shifting  each  detector  clockwise  (or  counterclockwise) 
by  one-quarter  of  a  detector  cell  width,  opposing  fan  beams  are  displaced 
laterally  by  one-half  of  a  cell  width  relative  to  each  other.  This 
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Figure  3.  Detector  Distance  Versus  Source  Distance  For 
System  II  and  III  Geometries. 


effectively  doubles  the  sampling  frequency  (see  sec.  3. 4. 1.1)  allowing  higher 
resolution  and  reduced  aliasing.^ 

Advantages  of  System  II  are  its  abl^ty  to  Incorporate  magnifications  In 
excess  of  2  more  easily  than  System  III  and  a  slightly  lower  shielding 
requirement  resulting  from  the  larger  effective  separation  between  tubes  and 
detectors  compared  to  System  III*  On  balance  however,  the  advantages  of 
System  III  significantly  outweigh  those  of  System  II. 

3.3  Electronic  Detector  Arrays’ 

The  basic  constraints  on  an  electronic  detector  array  for  the  BRL  Flash 
CT  System  are  similar  In  some  respects  and  different  In  others  to  those  of 
modern  medical  CT  scanners.  Ihe  x-ray  exposure  levels  per  view  and  the 
expected  dynamic  ranges  are  similar;  however  the  total  number  of  views  is  SO¬ 
SO  times  higher  for  medical  CT.  The  object  size  ranges  and  the  required  slice 
thicknesses  are  also  similar. 

On  the  other  hand,  the  required  spatial  and  contrast  resolutions  are 
lower  and  there  Is  no  limit  on  how  fast  the  data  must  be  read  out  even  though 
the  x-ray  pulse  time  Is  five  orders  of  magnitude  faster  than  medical  CT. 

Also,  there  can  be  several  centimetres  of  steel  present  In  objects  of  military 
Interest.  The  major  difference,  and  challenge  from  the  design  viewpoint.  Is 
the  much  higher  energy  x-ray  .beam  spectrum  In  the  flash  CT. 


The  higher  energy  x-ray  spectrum  has  several  ramifications.  First,  one 
needs  much  thicker  scintillator  crystals  to  absorb  the  Incident  radiation  and 
to  prevent  increased  noise  from  x-ray  photons  which  punch  through  the 
scintillator  and  then  Interact  In  rear-mounted  photodiodes.  Since  most  of  the 
lower  energy  photons  will  be  absorbed  via  photoelectric  Interactions  while 
most  of  the  higher  energy  x-rays  will  produce  Compton  scattering,  the  latter 
resulting  in  high-energy  outgoing  x-ray  photons,  these  scattered  x-rays  can 
then  Interact  In  adjacent  detectors  to  blur  the  resolution. 

As  a  result  of  all  these  considerations,  the  first  electronic  CT 
detectors  will  use  cesium  Iodide  scintillators  as  a  compromise  In  atomic 
number,  density  and  radiant  efficiency  (light  output).  Based  on  the  results 
obtained  with  this  scintillator,  future  electronic  detectors  may  use 
alternative  materials. 

The  anticipated  detector  design  parameters  are  listed  In  Table  2  (taken 
from  the  BRL  detector  specifications). 


^R.A.  Brooks,  G.H.  Glover,  A.J.  Talbert,  R.L.  Eisner,  F.A.  DlBlanca, 
"Aliasing-  a  Source  of  Streaks  In  Gomputed  Tomograms,"  J.  Oomput.  Assist. 
Tomography.  Vol.  3,  #4,  1979,  pp.  511-518. 
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Figure  4b.  Spatial  Resolution  Versus  Magnification  for  Various 

Detector  Slses 
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Figure  4e*  Schematic  Repreaentation  of  Zero  and  Quarter>-Offaet 
of  the  Detector  Oells  (Digitisation  Oells  for  Film)* 
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Figure  4d.  Effective  Nyquist  Resolution  Element  Versus 
Hagnification  for  Various  Detector  Sizes* 


Tftble  2.  Detector  System  Characteristics 


Detector  material 

width  by  height 
length 

radiation  efficiency 
Photodiode 


Collimator  plates 
dimension 


Cesium  Iodide  (Tl) 
0.302  by  2.5  cm 
2  cm 
>  12  Z 

CT-type,less  than  5Z 
variation 

In  radiant  sensitivity 
over  active  area 
Thngsten 

0.05  X  3.0  X  lA.S  cm 


Performance  parameters 
output  signal 

optical  crosstalk 
non-llnearlty 


average  over  5  x  10“^® 
coul/200  keV. 
x-ray  absorbed  In 
scintillator 

<  2  Z  between  adjacent 
cells 

<  2  Z 


3.4  System  Characteristics 


The  Image  quality  of  the  flash  CT  system  Is  characterized  by  three 
principal  parameters:  static  spatial  resolution,  dynamic  spatial  resolution 
and  contrast  resolution.  These  parameters  refer  respectively  to  the  system's 
ability  to  resolve  detailed  structure,  to  freeze  or  svold  blurring  In  highly 
dynamic  events  and  to  discriminate  between  two  adjacent  areas  which  differ 
slightly  In  x-ray  attenuation  coefficient. 


3.4.1  Static  Spatial  Resolution 

For  a  CT  scanner,  there  are  two  stages  at  which  the  spatial  resolution 
should  be  evaluated.  These  stages  are  after  the  view  projection  data  have 
been  accumulated  and  after  the  final  Image  has  been  reconstructed. 

3. 4. 1.1  Projection  Resolution 

The  spatial  resolution  of  an  Imaging  system  can  be  characterized  using 
the  modulation  transfer  function  where  o  Is  the  spatial  frequency  In 

units  of  length'^  and  represents  the  ratio  of  the  amplitude  of  the  output 
signal  to  that  of  a  sinusoidal  Input  signal  or  object  of  frequency  u.  The 
normalization  of  M(u)  is  usually  defined  such  thst  M(0)  ■  1. 

The  modulation  transfer  function  (also  referred  to  as  the  HTF)  equals  the 
modulus  of  the  Fourier  Transform  of  the  corresponding  "shape  fimctlon"  L(x) 
which  can  be  the  source  shape,  the  detector  response  shape,  a  combination  of 
the  two,  or  any  other  spatial  shape  function.  Thus, 

M(w)  -  |i*  L(x)  exp(-l(iK)  dx  j  .  (1) 

Since  the  combined  effect  of  two  shape  functions  Is  obtained  by 
convolutlng  the  Individual  shape  functions.  It  follows  from  the  above  equation 


that  the  combined  effect  of  two  MTF's  Is  obtained  by  multiplying  the 
Individual  MrF'8,(aee  Reference  6).  li>re  specifically,  the  total  MTF  of  the 
source  and  detector  Is  the  product  of  the  source  and  detector  HTF's  after  each 
Is  scaled  to  the  object  plane  using  the  system  magnification  m,  which  Is  the 
ratio  of  the  Image  size  to  object  size.  . 

Given  a  rectangular  shape  function  L(x)  with  half-width  h,  the  local 
(unsealed)  MTF  Is 

M(<i>)  -  sine  (  h  u  )  (2) 

where  sine  (  x  )  -  sin  (  x  )  /  x.  Plots  of  L(x)  and  M(u)  are  given  In  Figure 
4a.  The  total  MTF  at  the  object  plane  resulting  from  the  combination  of 
rectangular  source  and  detector  shapes  Is  then  (Reference  7): 

*tot^“^  -  slnc[(»-l)/m  hgQ  (d]  slnc[l/m  hj^^  u]  (3) 

where  the  source  and  detector  half-widths  are  h^^  and  h^.^  respectively.  The 
limiting  spatial  resolution  of  the  Imaging  system  Is  defined  here  as  the 
frequency  at  trtileh  the  system  MTF  drops  off  to  0.1,  l.e.  that  frequency,  m, 
where  the  output  amplitude  drops  off  to  lOZ  of  the  Input  amplitude.  Thus,  as 
shown  In  Figure  4a,  we  must  find  for  a  value  of  M  »  0.1.  The  limiting 
angular  spatial  frequency  can  now  be  used  to  determine  a  limiting  spatial 
resolution  p*  as 

p*  -  1/  «i^.  (4) 

The  equation  for  Is  transcendental  however,  l.e.,  at  cannot  be  found  In 
terms  of  so  the  limiting  value  of  at  (oF)  Is  found  by  Iteration  using  the 
Method  of  False  Position  and  the  result  Is  Indicated  in  Figure  4a.  p* 
corresponds  to  the  sampling  aperture  of  the  system  and  Is  shown  In  Figure  4b 
for  a  5  mm  source  width  and  several  values  of  detector  widths  and 
magnifications . 

Another  limiting  resolution  element  p»_  can  te  defined  corresponding  to 
the  Nyqulst  sampling  frequency  fj|y  of  the  detector^:  ' 

Pjly  -  1  /  f Ijy  -  2  h^et  ^5) 

If  objects  smaller  than  (l.e.,  frequencies  higher  than)  the  Nyqulst  frequency 
are  Imaged,  aliasing  artifacts  (streaks)  usually  result. 

In  using  P||y,  two  factors  must  be  considered.  First,  the  equation  must 
be  scaled  to  the'^object  plane.  Secondly,  the  detector  offset  from  the  system 


^T.H.  Newton,  D.G.  Potts,  Radiology  of  the  Skull  and  Brain.  TSchnlcal 
Aspects  of  Computed  Tomography.  The  C.V.  Mosby  Company,  St. Louis,  1981. 


'C.M.  Ooulam,  The  Physical  Basis  of  Madlcal  In 
New  York,  pp.  60,  97,  1981. 

^T.H.  Newton  and  D.G.  Potts,  (Ibid),  p.  3931. 
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centerline  mist  be  taken  Into  account  (see  Figure  4c) «  If  zero  offset  is 
used,  the  opposing  detector  pairs  will  align  and  the  Nyquist  frequency  will 
not  change.  If  "quarter  offset"  is  used,  opposing  fan  teams  will  be  offset  by 
one-half  the  sampling  distance  and  the  Nyquist  frequency  will  be  doubled 
(Reference  5).  The  detectors  in  System  III  opposite  one  another  are  nearly 
one-quarter  offset,  resulting  in  a  near-doubling  of  the  Nyquist  frequency. 


Combining  the  scaling  and 

offset  factors  gives: 

PNy  -  2  hjet  /  « 

(  0  offset  ) 

(6) 

PNy  ■  ^det  ^  ® 

(  1/4  offset  )  , 

(7) 

The  variation  of  pji^  with  detector  width  and  magnification  is  shown  in  Figure 
4d.  Inspecting  Figs.  4b  and  4d  leads  to  several  insights: 

1.  F6r  magnifications  only  slightly  less  than  two,  there  is  little 
gain  in  when  small  detector  elements  are  used. 

2.  When  the  detector  element  is  equal  to  (slightly  smaller  than)  the 
source  size,  the  optimum  magnification  is  2  (between  about  1.6  and  2);  but 
when  a  very  high  resolution  detector  such  as  screen  film  is  used,  the 
resolution  continues  to  improve  down  to  very  low  magnifications.  These 
conclusions  are  straightforward  results  of  Equation  3. 

3.  A  quarter-offset  alignment  should  theoretically  help  to  reduce 
aliasing  and  thereby  allow  higher  spatial  resolution.  In  practice,  however, 
the  deleterious  effects  from  the  severely  limited  number  of  views  (view 
aliasing)  may  eliminate  most  of  this  benefit. 

4.  As  magnification  increases,  puy  continually  improves  since  the 
sampling  distance  at  the  object  plane  continually  decreases. 


As  a  check  on  the  equations,  the  theoretical  values  for  limiting  spatial 
resolution  were  calculated  and  compared  to  published  values  (Reference  9)  for 
two  commercial  CT  scanners  and  found  to  agree  within  15-20  Z.  These  machines 
are  not  view  limited  as  reconstructed  images  are  formed  from  hundreds  of- 
projections  rather  than  the  21  used  in  the  BRL  system.  However  one  expects 
the  CT  resolution  of  the  BRL  design  to  be  degraded  from  the  theoretical 
resolution  by  more  than  this  amount.  On  the  other  hand,  the  projection 
resolution  of  the  BRL  machine  should  be  close  to  the  theoretical  estimates. 

3. 4. 1.2  Image  Resolution 

The  image  resolution  is  more  difficult  to  obtain  than  the  projection 
resolution.  The  image  resolution  is  always  degraded  from  the  projection 
resolution  by  an  amount  which  depends  on  the  number  of  views  available  and  on 
the  algorithm  used  for  reconstruction.  Thus,  the  final  resolution  can  be 


^Anonymous,  "CT/T  8800  Technology  Brochure,  General  Electric  Medical  Systems," 
1979;  "CT  9800  Evaluation  Criteria,"  General  Electric  Medical  Systems,  1981. 
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Figure  6.  Schematic  of  Water  Phantom  Used  for  Projection  and  CT 
Contrast  Resolution  Determination* 
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determined  only  with  the  actual  reconstruction  algorithm  operating  on  real  or 
accurately-simulated  data  and  Is  not  given  here. 

3.4.2  Dynamic  Spatial  Resolution 

The  dynamic  spatial  resolution  refers  to  Image  blurring  from  object 
motion  perpendicular  to  the  beam  direction  during  the  x-ray  pulse  and  can  be 
characterized  by  an  appropriate  MTP.  The  effect  of  such  object  motion  can  be 
expressed  as  a  convolution  of  the  system's  static  line  spread  function  (LSF) 
with  a  dynamic  LSF.  The  shapb  of  the  dynamic  LSF  will  be  the  same  as  that  of 
the  x-ray  intenslty-versus-tlme  curve  but  having  a  width  equal  to  the  object 
motion  during  the  x-ray  pulse.  It  Is  convenient  to  determine  the  MTF 
corresponding  to  the  motion  blurring  so  that  the  total  system  MTF  Including 
motion  can  be  expressed  as  the  product  of  the  static  and  dynamic  NTF's. 

Since  the  x-ray  pulse  time  Is  approximately  25  nsec,  the  corresponding 
dynamic  MTF  for  projection  radiography  Is  virtually  perfect  (flat  with  a  value 
of  one)  for  all  spatial  frequencies  of  Importance  in  ballistics  studies. 
However,  for  CT  operation  one  cannot  assume  that  all  21  sources  will  fire 
simultaneously.  Calling  the  time  over  i^ich  all  tubes  discharge  t  and  the 
velocity  of  the  moving  structures  v.  the  dynamic  MTF  is  given  by: 

M(u)d  -  slnc(  vt(V2).  (8) 

Assuming  that  t  >  1  mlcrosec  and  v  -  1  and  10  Im/sec  yields  a  dynamic  MTF 

as  shown  In  Figure  5.  At  velocities  of  up  to  1  km/ sec  there  is  little  loss  of 
resolution  for  objects  larger  than  2  mm. 

3.4.3  Contrast  Resolution 

As  was  done  for  spatial  resolution,  one  can  also  consider  the  contrast 
resolution  In  each  projection  as  well  as  that  of  the  final  reconstructed  CT 
Image.  Both  are  treated  In  this  section. 

3. 4. 3.1  Projection  Contrast  Resolution 

To  calculate  the  contrast  resolution  In  each  projection,  we  consider  a 
uniform  water  phantom  of  thickness  L  containing  a  small  object  having  the  same 
mass  attenuation  coefficient  as  water  but  a  different  electron  density.  The 
object  Is  taken  to  be  a  cylinder  with  a  diameter  of  1  cm.  1cm  In  length  and 
with  a  density  slightly  different  from  water.  This  configuration  Is  shown  In 
Figure  6 . 

The  projected  density  contrast  C  Is  defined  as  the  fractional  difference 
In  the  density  of  the  small  cylinder  D  and  the  background  D' . 

C  -  (D  -  D')  /  (.5  X  (D  +  D')).  (9) 

The  x-ray  energy  spectrum  was  calculated  on  an  IBH-PC  computer  with  the 
program  XRSPBC  which  determines  the  bremsstrahlung  and  characteristic  spectra 
for  an  arbitrary  voltage  and  current  waveform.  The  waveforms  for  the  1  Wf 
flash  x-ray  source  used  at  BRL  were  obtained  by  scaling  the  data  In  Figure 


2.32  of  Reference  10  to  a  voltage  of  1  MV.  The  program  also  calculates  the 
exposure  for  a  given  electron  charge  transported  In  the  x-ray  tube  (called  mAs 
or  current-time  product  when  the  current  and  voltage  are  constant) . 


The  detector  was  modeled  as  a  cesium  Iodide  crystal  array  with  a  2  cm 
thickness.  ■  The  signal  was  taken  as  the  total  energy  absorbed  In  the 
scintillator  excluding  reabsorption  of  scattered  x-rays.  The  Image  noise  was 
considered  to  arise  exclusively  from  the  statistical  fluctuations  In  the 
absorbed  energy.  Furthermore,  the  x-ray  flux  was  scaled  to  agree  with 
experimental  measurements  taken  with  a  4  mm  focal  spot  flash  x-ray  tube  and 
then  adjusted  for  the  2.1  m  source-to-center  distance. 

The  contrast  resolution  0*  Is  defined  as  the  threshold  density  contrast 
which  produces  a  slgnal-to-nolse  ratio  of  8.  This  value  Is  based  on  the 
typical  results  of  observer  stiidles  In  medical  x-ray  Imaging  and  represents 
the  smallest  perceivable  slgnal-to-nolse  ratio.  As  an  example.  If  the  density 
contrast  of  two  objects  can  be  resolved  by  an  observer  then  the  difference  In 
the  grey  scale  value  of  the  two  objects  divided  by  the  average  noise  In  the 
Images  will  be  8,  or  larger.  For  a  single  X!-ray  pulse,  the  resulting 
projection  contrast  resolution  Is  given  In  Table  3. 

Table  3.  Calculated  projection  contrast  resolution  versus  phantom  diameter. 


L(cm)  C*(Z) 

01  0.7 

10  1.4 

20  2.7 

30  5.3 

3.4. 3. 2  Image  Contrast  Resolution 

The  contrast  resolution  of  the  final  CT  Image  may  be  estimated  by  two 
different  methods  given  In  this  section.  The  first  method  uses  the  projection 
contrast  resolution  calculated  In  the  previous  section,  while  the  second 
method  Involves  scaling  published  data  from  a  medical  CT  scanner.^*  Slhce  the 
latter  data  were  obtained  with  a  16.5  cm  diameter  water  phantom  and  a  1  cm 
slice  thickness,  the  same  values  are  used  here. 

Assuming  that  the  CT  contrast  resolution  scales  with  the  total  Image 
noise  (l.e.,  with  the  Inverse  square  root  of  the  total  number  of  x-ray  photons 
detected  In  all  views)  and  performing  a  log-log  Interpolation  of  the  data  In 
Table  3  to  an  L  value  of  16.5  cm  giving  O*'  -  2.3.  Since  In  the  BRL  System  III 
there  are  21  projections. 


^^F.  Jamet,  G.  Ihomer,  Flash  Radiography,  Elsevier  Scientific  Publishing 
Co.,  Amsterdam,  1976. 

Cohen,  F.A.  DlBlanca,  "The  Use  of  Oontrast-Detall-Dose  Evaluation  of 
Image  Quality  In  a  Computed  Tomographic  Scanner,"  J.  Oomput.  Assist. 
Tomography.  Vol.3,  #2,  pp.  189-195,  1979. 
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CcT*-  2.3/(21)^/2j;  .  q.S  X  (method  1) 

For  the  second  method,  which  also  assumes  that  the  contrast  resolution 
varies  Inversely  with  the  square  root  of  the  total  number  of  x-ray  photons 
detected,  we  begin  by  estimating  the  total  radiation  exposure  delivered  In  a 
21  view  CT  scan.  Summing  the  radiation  exposure  along  the  lines  from  each 
source  position  to  a  point  on  the  phantom  surface  and  attenuating  by  the 
appropriate  chord  lengths  through  the  phantom  gives  an  estimated  total 
exposure  of  233  mR.  Then  the  CT  contrast  resolution  can  be  estimated  from 

CcT*  =■  k  (  R  h  (11) 

where  R  Is  the  exposure,  h  Is  the  slice  thickness  and  k  Is  the  unknown 
proportionality  constant.  Using  the  values  from  Reference  11,  namely  C*  > 

0.16  Z,  R  -  5000  mR  and  h  >  1  cm,  we  determine  k  and  then  obtain  for  the  BRL 
flash  CT  system 

CcT*  ~  0.8  %  (method  2).  (12) 

Considering  the  widely  differing  approaches  and  assumptions  (except  for 
the  basic  assumption  related  to  the  total  it-ray  detection) ,  many  of  which  are 
only  approximations,  used  to  obtain  these  estimates  of  CT  contrast  resolution, 
the  agreement  between  the  two  methods  Is  rather  remarkable.  These  estimates 
lead  us  believe  that  the  BRL  flash  CT  scanner  will  have  an  actual  contrast 
resolution  near  IZ  for  1  cm  diameter  pins  In  a  16.5  cm  diameter  water  phantom 
with  a  1  cm  slice  thickness.  The  above  analysis  does  not  evaluate  deleterious 
effects  of  the  artifacts  expected  in  the  BRL  system,  for  example  artifacts  due 
to  limited  number  of  views. 


4.  EXPERIMENTAL  STUDIES 


4.1  Background 

Previous  studies^  on  the  feasibility  of  a  ballistic  tomography  system 
were  concerned  with  the  quality  of  a  cross-sectional  Image  produced  with  a 
limited  number  of  views  or  projections.  Since  between  15  and  25  were 
considered  to  be  the  maximum  practical  number  of  views,  the  impact  of  this 
limitation  on  the  Image  was  considered  to  be  the  most  Important  aspect  of  the 
study.  Consequently,  the  rest  of  the  experimental  set-up  was  Idealized  so  as 
to  minimize  other  problems.  An  essentially  monochromatic  source  was  used 
along  with  a  well  collimated  source  and  detector  so  that  scattering  was 
minimized.  A  parallel-beam  geometry  was  also  used  which  simplified  the 
algorithm.  As  Is  seen  In  Figure  1,  the  real  system  Is  somewhat  more 
complicated.  The  flash  x-ray  source  required  to  record  the  event  In  the 
microsecond  time  frame  Is  not  monochromatic  but  has  a  bremsstrahlung  spectrum 
starting  at  low  energies  and  extending  up  to  the  maximum  discharge  voltage 
which  can  be  as  high  as  1  MV.  Since  the  source  Is  a  single  pulse  system  the 
dose  output  Is  limited  and  cannot  be  Integrated  in  time.  The  colllmatlon  will 
not  be  Ideal  with  a  fan  beam  geometry  and  there  will  be  a  scattering 
contribution  to  each  projection.  The  focal  spot  size  for  some  flash  x-ray 
systems  can  be  as  large  as  5  mm.  This  will  lead  to  a  blurring  of  the  Image  on 
each  projection.  The  fan  beam  geometry  will  require  a  modified  algorithm. 
Moreover,  the  radiation  field  produced  by  the  flash  x-ray  source  must  be 
uniform  over  the  object  or,  at  least,  reproducible  otherwise  accurate 
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attenuation  measurements  will  not  be  possible.  The  first  detector  that  will 
be  used  will  be  medical  x-ray  film  and  Intensifying  screens,  rather  than  a 
scintillator  and  photodiode  detector.  This  avoids  the  necessity  of  having 
over  2000  Individual  scintillators  and  associated  data  acquisition 
equipment.  It  also  allows  for  the  possibility  of  acquiring  conventional  two- 
dimensional  radiographs  of  ballistic  events.  Also  screen-film  systems  have 
high  spatial  resolution.  However,  film  has  a  relatively  limited  dynamic  range 
with  respect  to  Intensity  and  also  does  not  have  a  wide  range  of  linearity. 

The  quantum  detection  efficiency  Is  significantly  less  than  the  single 
crystal,  photodiode  array,  leading  to  Increased  Image  noise  and  reduced 
contrast  resolution.  All  of  the  above  differences  and  uncertainties  must  be 
examined  to  determine  the  Impact  on  the  reconstructed  Image. 

4.2  Single  Channel  Experiments 

The  source  used  In  acquiring  the  data  was  a  Hewlett-Packard  1  flash  x- 
ray  system.  The  propellant  mock-up  was  placed  Inside  a  fiberglass  chamber 
which.  In  turn,  was  located  Inside  a  lead  container.  Radiographs  were  taken 
through  a  horizontal  silt  In  the  lead  container.  A  photograph  of  the  side 
view  of  the  propellant  mock-up  Is  shown  In  Figure  7a.  A  typical  radiograph  Is 
seen  In  Figure  7b.  An  Intensity  profile  of  the  radiograph  Is  shown  In  Figure 
7c. 

This  x-ray  source  delivers  a  nominal  dose  of  65  mR  at  a  distance  of  1 
metre  and  a  nominal  focal  spot  size  of  5  mm  with  a  pulse  duration  of  25 
nanoseconds.  The  focal  spot  size  was  measured  by  taking  a  radiograph  of  a 
tungsten  plate  and  measuring  the  edge  gradient  as  recorded  on  film.  The  dose 
was  measured  with  a  Radcal  Model  lOlSC  dosimeter  which  has  an  accuracy  of 
1%.  The  reproducibility  of  the  source  was  determined  by  measuring  the 
deviation  of  the  dose  over  twenty  one  pulses.  The  standard  deviation  was 
measured  to  be  2X.  The  x-ray  field  uniformity  was  measured  over  a  beam 
(plane)  angle  of  6*  which  covers  almost  any  object  of  current  Interest.  A 
medical  film-screen  cassette  was  placed  at  a  distance  of  2  metres  from  the 
source  and  was  subjected  to  one  pulse.  The  film  was  carefully  developed  In  an 
automatic  processor  so  as  to  avoid  developing  artifacts.  The  optical  density 
of  the  film  was  measured  with  an  X-Rlte  Model  301  densitometer  with  a 
resolution  of  ±  0.01  D.  The  optical  density  measured  was  In  the  linear  region 
of  the  film  (1.52  D) .  The  variability  of  the  optical  density  was  less  than  ± 
0.02  D  (or  ±  2%  source  Intensity)  over  the  area  of  the  film  corresponding  to 
the  6"  angle.  The  source  was  heavily  shielded  with  lead  plates.  It  was  found 
that  room  scattering  was  a  significant  problem,  giving  unwanted  scattered 
radiation  at  the  detector.  The  shielding  was  required  to  reduce  this 
problem.  The  average  energy  of  the  source  was  found  by  measuring  the 
attenuation  through  several  thicknesses  of  Iron.  The  slope  of  the  plot  of 
ln(l)  versus  thickness  yields  an  effective  attenuation  coefficient  from  which 
the  average  energy  Is  obtained  from  cross  section  tables.  Measurements  made 
with  both  lead  and  Iron  yielded  an  average  energy  of  250  keV.  The  beam  Is 
hardened  after  passing  through  0.64  cm  of  steel  and  the  emerging  radiation  has 
an  average  energy  of  600  keV. 

The  model  seen  In  Figure  7a  was  mounted  on  an  Ealing  Model  35-2500  rotary 
stage.  Eighteen  projections  were  used  In  the  reconstruction  and  the  rotary 
stage  was  used  to  accurately  position  the  object.  The  stage  has  an  accuracy 
of  better  than  one  minute  of  arc  and  a  deviation  from  horizontal  plane  of  less 


than  one  arc  minute.  It  can  support  a  weight  of  up  to  25  kg.  For  the 
projection  data  to  be  accurate,  the  axis  of  rotation  must  be  perpendicular  to 
x-ray  radiation  plane.  To  assure  this,  four  positioning  arrows  were  attached 
vertically  to  the  rotation  table.  Radiographs  were  recorded  of  the  four 
arrows  at  0**,  90°  and  180°  Intervals.  The  ends  of  the  arrows  should  form  a 
straight  line  on  the  radiograph.  The  rotating  table  was  adjusted  until  this 
was  achieved.  An  example  of  the  radiograph  Is  shown  In  Figure  7b.  By  having 
a  fixed  object  In  the  radiograph,  the  axis  of  rotation  can  be  determined  from 
the  0°  and  180°  positions. 

The  detector  used  for  these  tests  was  a  cassette  of  3M  Trlmax  12  rare 
earth  screens  and  XM  film.  This  combination  has  a  very  high  sensitivity  with 
a  slight  loss  In  resolution  over  more  conventional  combinations.  However,  as 
was  determined  In  Section  3. 4. 1.1,  spatial  resolution  In  the  projection  was 
not  a  major  factor  In  the  reconstructed  Image  but  slgnal-to-nolse  ratio  was 
found  to  be  very  Important.  Additionally,  the  film  response  Is  more  linear 
when  the  optical  density  Is  In  the  0.7  to  2.5  region.  This  requires 
significant  sensitivity  for  the  Intensifying  screens.  The  question  of 
linearity  of  the  detector  system  remains.  The  system  was  calibrated  using  six 
plastic  rods  of  different  length,  each  one  shielded  with  a  lead  sleeve  so  that 
scattering  from  one  would  not  contribute  a  signal  to  the  other.  These  were 
recorded  on  the  same  film  with  the  projection  data.  Examination  of  the 
calibration  data  of  optical  density  versus  length  Indicated  that  the  system 
was  linear  over  the  optical  densities  of  Interest  hence  recallbratlon  of  the 
data  was  not  required. 

4.3  The  Data  Set 

Data  were  taken  with  a  source-to-object  distance  of  3.31  metres,  and  with 
an  object-to-detector  distance  of  0.223  metres.  This  configuration  minimized 
the  problem  associated  with  the  large  (5  mm)  focal  spot  size.  The  vertical 
height  of  the  Image  represents  the  collimatlon  window  which  was  25  mm. 

Optical  density  profiles  of  these  data  were  acquired  In  two  ways.  First  an 
Optronics  microdensitometer  was  used  employing  a  0.2-mm  sampling  orifice  to 
digitize  the  radiograph  numerically.  Six  horizontal  lines  were  averaged  to 
produce  a  density  profile  over  a  height  of  1 .2  mm.  These  data  were  used  In  the 
reconstruction  algorithm.  An  example  of  such  a  profile  Is  given  In  Figure 
7c.  An  Eyecom  vldlcon  camera  and  a  Comtal  Vision  One/20  Image  analyzer  were 
also  used  to  produce  a  series  of  profiles.  Data  were  also  taken  with  other 
object  and  detector  distances,  but  they  have  not  yet  been  used  to  reconstruct 
Images . 


5 .  RECONSTRUCTED  IMAGES 


The  Algorithm 


The  number  of  views  available  from  the  BRL  facility  Is  only  one  tenth 
that  of  a  conventional  CT  system.  Consequently,  considerable  emphasis  was 
placed  during  the  design  phase  of  the  facility  In  arriving  at  the  selection  of 
a  reconstruction  algorithm  which  produces  the  least  artifact- prone  Image. 
Computer  experiments^  showed  that  maximum  entropy  reconstructions  were 
superior  to  conventional  algorithms  such  as  filtered  back  projection  as  far  as 


Image  quality  waa  concerned.  The  higher  quality  though  Is  bought  at  the 
expense  of  conslderably^lncreased  computer  running  time  requirements. 

Maximum  entropy  Is  based  on  several  premises.  This  Includes  the  notion 
that  the  Image  with  the  minimum  Information  Is  the  one  represented  by  maximum 
entropy.  Thus  the  object  which  transmits  the  minimum  cunount  of  Information 
about  Itself  has  to  be  the  one  selected.  Further,  maximum  entropy  Is  valid 
without  modification  only  when  the  object  consists  of  isolated  points.  The 
maximum  entropy  Image  contains  only  the  structure  embodied  In  the  data 
consequently  leading  to  a  minimization  of  the  number  of  artifacts  present  In 
the  Image. 

An  already  existing  computer  code,  MEMT,  was  chosen  for  the  data 
reconstruction.  A  photograph  of  a  mock-up  along  with  the  reconstructed  Image 
using  the  fan  beam  MENT  algorithm  are  seen  In  Figure  8.  The  reconstruction 
gives  a  reasonable  representation  of  the  mock-up.  MENT^^  Is  based  on  the 
concept  of  maximum  entropy  discussed  at  length  In  Reference  13,  and  Involves 
the  following  steps:  transmission  data  from  the  experiments  are  available  In 
the  form 

®jm  “  //f(x,y)dxdy  (13) 

where  are  the  projection  data,  l.e.,  Ind/I^),  the  ratio  of  the  measured 
to  the  Incident  radiation,  and  f(x,y)  the  attenuation  In  the  region  of  the 
object  under  study.  Since  many  functions  satisfy  this  equation,  a  method  was 
devised  to  select  the  function  which  came  closest  to  picking  out  the 
attenuation  observed.  For  this  a  second  function,  called  the  entropy  Is 
defined  where 

C  -  //f(x,y)  ln(f(x,y)  A  )  dxdy  (14) 

and  a  constrained  optimization  problem  Is  solved.  That  Is,  one  must  find  f 
which  maximizes  (  subject  to  the  constraints,  l.e.  the  measured  values.  It  Is 
known  that  for  a  maximum  to  exist,  Lagrange  multipliers  must  exist  such  that 
f(x,y)  Is  a  stationary  point  of  the  Lagrangian.  For  this  to  be  true,  the 
functional  derivative  of  the  Lagrangian  must  be  zero,  leading  to  a  system  of 
nonlinear  equations  which  are  solved  Iteratively  by  the  Gauss-Seldel 
technique.  The  solution  converges  to  the  image  with  the  least  amount  of 
Information,  l.e.,  the  one  which  has  the  fewest  number  of  artifacts. 

5.2  Results 

In  Figure  9,  a  reconstructed  image  of  a  mathematical  mock-up  consisting 
of  a  solid  lucite  rod  In  the  center,  an  aluminum  cylinder  filled  with  water  on 
the  left,  and  a  lucite  tube  filled  with  water  on  the  right  are  shown.  The 


^^6.  Hlnerbo,  ”MENT:  A  Maximum  Entropy  Algorithm  for  Reconstructing  a  Source 
from  Projection  Data,"  Computer  Graphics  and  Image  Processing.  Vol.  10, 
pp.  46-68,  1979. 

Skilling,  "The  Maximum  Entropy  Method,”  Nature ,  Vol.  309,  pp.  748-9, 
1984. 


density  difference  between  the  Incite  and  the  water  is  apparent.  The  spatial 
resolution  capability  of  the  system  Is  Illustrated  In  Figure  10  in  a  mock-up 
consisting  of  an  aluminum  tube  containing  Incite  rods  of  5*  3»  and  1-mm  In 
diameter.  The  5  mm,  and  even  the  3  mm  and  1  mm  rods  are  clearly  visible. 
Reconstruction  times  are  typlbally  of  the  order  of  8  minutes  on  a  Cyber  76 
computer. 

6.  CONCLUSIONS 

.  The  first  microsecond  tomographic  facility  under  construction  at  the 
Ballistic  Research  Laboratory,  offers  a  four  orders  of  magnitude  Improvement 
In  the  time  resolution  capability  over  the  current  generation  of  scanners 
(Figure  11).  With  a  spatial  resolution  of  2  to  4  millimetres  on  an  object  of 
200  mm,  a  density  resolution  of  better  than  IZ  and  time  resolution  In  the 
microsecond  range,  multidimensional  density  profiles  of  transient  events  of 
ballistic  Interest  will  be  amenable  to  study,  nonlntruslvely  for  the  first 
time.  We  believe  that  despite  the  high  source  voltages,  adequate  shielding 
can  be  provided  against  Compton  scattering.  The  facility  can  also  be  used  for 
limited-view  tomographic  clneswtography.  These  conclusions  did  not  consider 
view  Halted  artifacts  and  It  Is  recognised  that  system  performance  may  be 
degraded  under  actual  firing  condltlona.  It  Is  also  Interesting  to  note  that 
the  time  resolution  was  obtained  without  any  spatial  resolution  loss. 
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Figure  11.  Scan  Times  of  the  Current  Generation  of  Scanners 
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APPENDIX  A.  EQUATIONS  FOR  SYSTEM  GEOMETRIES 

The  equations  for  the  most  compact  geometries  for  Systems  II  and  III  are 
derived  here.  We  begin  with  System  III  because  the  mathematical  treatment  Is 
simpler.  System  II  requires  an  Iterative  solution;  It  cannot  be  solved  In 
closed  form.' 


A.l 


Figure  12  shows  the  basld  geometry  of  System  III  In  Its  most  compact  form 
In  uhlch  each  of  two  nearly  opposite  x-ray  beams  Just  clear  each  other's 
detector.  The  major  variables  are: 


dg  ■  source-to-centefr  distance 
dQ*detector-to-center  distance 

d  «  2r  "  scan  field  diameter  (largest  object) 

^  -  half  angle  of  fan  beam 

a  "  half  angle  between  adjacent  sources  • 

By  straightforward  trigonometry  one  has: 

tan  0  -  X  /  [  (dg  +  djj)  sin  ^  ]  (Al) 

and  X  ■¥  X  !  cos  6  -f  dg  *  (djj  ■¥  dg)  cos  #  .  (A2) 

Substituting  6"x/2-  ^  -  a,  combining  the  above  two  equations  and  solving 
for  dQ  In  terms  of  dg,  d  and  Ng  gives: 

dp  -  I  d/2  +  dg  sin  (  ♦  +  a  )  1  /  sin  o  -  dg  (A3), 

(Here,  a  ■>  «  /  Ng  and  t  **  arcsln  (  0.5  d  /  dg  ).  Further,  we  note  that  the 
number  of  sources  Ng  must  be  odd  for  a  maxlmally-compact  System  III.  This  Is 
not  the  case  for  System  II. 


System  II  Geometry 


The  most  compact  geometry  for  System  II  Is  given  In  Figure  13.  Here, 
adjacent  detectors  must  Just  clear  each  other.  The  geometrical  variables  are 
defined  similarly  to  those  for  System  III  except  that  a  now  represents  the 
full-angle  between  the  adjacent  sources. 


System  II  geometrical  equations  are  Identical  to  Equations  Al  a'nd  A2  but 
now  0-Tr/2+#-a.  Solving  for  dp  gives: 

d„  -  ld/2  +  dg  sln(a  -  ♦  )I  /  sln(o-  2*  )-dg  .  (AA) 

Here,  ^  Is  defined  as  It  was  for  System  III  but  a  complication  arises  In 
the  definition  of  a.  For  System  II,  the  angle  between  sources  Is  limited  by 
the  necessity  that  the  beams  at  either  end  of  the  source  arc  clear  each 
other's  detectors.  This  Is  precisely  the  constraint  arising  In  System  III. 
Therefore: 


Ojj  -  (  v  -  ojjj  )  /  (  Ng  -  1  )  (A5) 
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where  Is  the  solution  to  Equation  A3;  but  is  as  yet  undetermined! 

Solving  Equation  A3  for  gives: 

slno^T  -  lAd/2+dgSln^A^  - 

d2/4+dg2 

where  A^  dp  +  dg  (  1  -  cos  4)*  By  inspecting  trial  solutions  of  the 
quadratic  Equation  A6,  one  finds  that  the  positive  square  root  is  correct. 

Now  we  can  find  dQ  for  System  II  by  the  following  iterative  procedure: 

1.  seto^^j  to  zero; 

2.  solve  Equation  AS  for 

3.  solve  Equation  AA  for  dp  -  if  converged,  go  to  5; 

4.  solve  Equation  A6  for  using  dp  then  loop  to  2; 

5.  finished  -  record  final  value  for  dp. 

Solutions  for  both  system  geometries  are  Implemented  on  an  IM-PC  computer  in 
a  BASIC  program  called  BRLGEOM  which  executes  in  a  few  seconds. 
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